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Abstract 

A logical explanation as to why the choice of 

1 t? 
1 f 
10 

(the Pian-Sumihara basis) as a linear basis to approximate stress leads to greater 
efficiency in enhanced strain problems, is presented. An Airy stress function and the 
consequent selective simplification resulting from the differentiation of an implied, 
single, parent approximating polynomial, are the essence of this argument. 
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1 Introduction 



Pian and Sumihara first identified the basis 

" 1 n " 

1 £ 
_ 1 _ 

as the most efficient linear basis for approximating stress in enhanced strain problems. 
This observation they made more rigorous by way of a Wilson element (a perturbation 
of sorts). 
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This paper presents a logical mathematical argument for making the same choice of basis, 
albeit with the wisdom of hindsight. It attributes the greater efficiency of the basis to 
properties inherent in the mathematics of the problem. The components of the stress 
tensor are recognised to be related by way of an Airy stress function and it is in this way 
that a fundamentally more correct representation of the full linear basis is arrived at. By 
further desiring the advantages of a two field problem, the most efficient, linear basis is 
obtained. 



2 An Airy Stress Function 

The Airy stress function is a potential of sorts. Interpreting stresses to be the various sec- 
ond derivatives of a single polynomial leads to selective simplification and interdependence 
between the resulting linear approximations. This simplification and the interdependence 
are not obvious in a more superficial treatment. 

n . da u da 12 da 21 da 22 
div cr = =^ — h -7T- = and 1" — = 

ox oy ox oy 

This is recogniseable as 

curl(-<7i2, (Tn, 0) = and curl(<7 2 2, —021,0) = 0. 

This, in turn, implies that (— 012,011, 0) and (0-22,-02150) may be interpretted as Va 
and V/9 respectively, without any inconsistancy in the 

curlV( • ) = 

identity. 

By symmetry of cr, 

da 8(3 

012 = 021 T\ = 

ox oy 

and for a two dimensional problem of the type under consideration this once again implies 

curl(/?,a,0) = 0. 

(/?, a, 0) may therefore be interpretted as V$ without any inconsistancy in the 

curlV( • ) = 

identity. 

In summary, with an equation 

div cr = 

governing the motion, in the two-dimensional case, the components of the stress may be 
derived from an Airy stress function as follows 
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022 

where $ is the Airy stress function. 



d 2 $ 

dx 2 ' 
<9 2 $ 

dxdy ' 



2.1 Finite Element Approximation 

Due to approximation, 

div cr = 

and not the constitutive 

div a = f 

are really the equations being solved (Reddy 0). 
Defining a function 

on each element Q e , 

<9 2 $ 



022 



<9x 2 

d ( d(j) d^ d<p drj \ d£ t d ( dtp d£ dcj) drj \ drj 



d£ \ dt; dx\ drj dx\ J dx\ drj \ d£ dx\ dr] dx\ ) dx\ 
I d 2 4> d£ d 2 4> drj \ <9£ ( d 2 (j) <9£ d 2 <f> dr] \ dr] 
\d^ 2 dxi dr/d!; dx\ J dx\ \ dr/d^ dx\ drf dx± J dx\ 



Assumption 

The individual elements, Q e , are usually mapped to the master element, Q, with ~ 
J^- ps on average, J^- » a\ and -j^- ps a,2, a>i and a>2 some constants, on average. 
(Alternatively it can be argued that there will be no loss of generality or weakening of 
the argument if a rectangular mesh is considered. Not allowing this simplification leads 
to an extremely messy argument, a chapters long exercise in differentiation.) This implies 



022 = «l 



2 d\ 



Similarly, 

011 



di 2 

d 2 6 



2 

2 Q^2 

d 

0~12 = CT 2 i = CLlCl2 



d^drj 
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3 The Relationship Implicit in the Linear Approxi- 
mation 

Since linear approximations of an are to be considered, 

&1 + »24 + M 



where b%, b 2 and 63 are the relevant combining constants. This means 



0(0 ??) = J i h + b 2 £ + hi] drjdrj 



= c 1 + c 3V +h 1 r ] 2 + ^ V 2 + ~h v 3 + vh(0+f2(0 (i) 

in which the exact form of 77/1(77) + ^(77) remains to be determined. Similarly, approxi- 
mating cr 2 2 as some multiple of 64 + + b§r] implies this very same polynomial function 

0(0 = J J ^4 + ^£ + ^6^7 (by Airy stress function) 

= c 1 + c 2 ^ + h 4 e + h 5 e + lb 6 ev + ^9i(v)+92(v), (2) 

ZD/ 

in which the exact form of (72(77) is determined by equation (|I|). This equation in turn 
specifies /2(C) m equation flU). Approximating cri 2 = cr 2 i in it's turn as as 67 + b${; + 6977 
implies the polynomial function 

el r l 



</>(£, 77) = J_ J br + bst + bgvdtdri 



= c 1 + b 7 ^ V +^b 8 ev + lb 9 ^ 2 + f 2 (0 + 92(v) (3) 

where /2(C) an d g 2 (v) have already been determined by equations @ and (]!]) respectively. 
This last expression for (/>(£, 77) also specifies the, until now undetermined, 77/1(0 and 
£(71(77) in equations ([!]) and (|2|). In summary, collecting equations ([!]), (0) and (|3]) together 
leads to the specification of an implied, single parent approximating polynomial 

0(0 ?7) = Ci + C 2 £ + C377 + C 4 + C 5 ^77 + C 6 77 2 + C 7 + C 8 0t? + Cg^T/ 2 + CioT/ 3 . 

Having established both the existance and nature of the relationship between the con- 
stants in what were apparently seperate linear approximations, 



2 



di] 2 
d 2 <p 
d^drj 



2c 4 + 6c 7 £ + 2c 8 77 
2c 6 + 2c 9 £ + 6ci 77 
c 5 + 2c 8 £ + 2c 9 77 



can now be written where the Cj's (i = 4, ■ • • 10) are constants related to the finite element 
solution of the problem in question. 
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Conclusion 

The Airy stress function therefore reveals how a linear approximation of the components 
of cr on each element really amounts to 
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instead of the superficially more obvious 
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4 Eliminating the Last Two Columns 

The rank of the matrix in equation (|4]) indicates that there are still two extra columns. 
The equation in which cr is used is a three-field problem, in which the strain, 7, only 
occurs once in a term cr-j. Choosing cr correctly would reduce the problem to a two-field 
problem since 



J cr ■ -fold = 



is required in accordance with Reddy H. In other words 



cr • 7 



1 r/ f 
1 f r] 
1 i] i 



£000 
77 
£ 77 



must always be zero. This is only certaint if the sixth and seventh columns of the stress 
basis are omitted. 



5 Conclusion 



An Airy stress function and consequent simplification resulting from the differentiation 
of an implied, single, parent, approximating polynomial are able to provide a logical 
explanation as to why the choice of 

" 1 77 " 
1 i 
10 



(the Pian-Sumihara basis) as a linear basis to approximate stress leads to greater effi- 
ciency in enhanced strain problems. 
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